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Abstract 

The influence of nonvanisliing quark masses on tlie total cross section 
in electron positron collisions and on the Z decay rate is calculated. 
The corrections are expanded in m^/s and as- Methods similar to 
those applied for the quadratic mass terms allow to derive the cor- 
rections of order a^m^/s^ and a^m^/s^. Coefficients which depend 
logarithmically on m^/ s and which cannot be absorbed in a running 
quark mass arise in order a^. The implications of these results on elec- 
tron positron annihilation cross sections at LEP and at lower energies 
in particular between the charm and the bottom threshold and for 
energies several GeV above the bb threshold are discussed. 
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1 Introduction 



The precise measurement of the cross section for hadron production in electron 
positron annihilation leads to an accurate determination of the strong coupling 
constant a^. The method is free from a variety of theoretical uncertainties and 
ambiguities which beset many other observables of (still) superior statistical ac- 
curacy. The extraction of from the ratio R between the rates for hadron and 
lepton pair production is based on the calculation of perturbative QCD correc- 
tions. These are available up to third order in in the massless limit |1|, 0] . Mass 
corrections were calculated for small m? / s up to order for vector and up to 
order a1 for axial vector current induced reactions ^. In these papers it was 
furthermore demonstrated that all large logarithms which arise in the expansion 
coefficients can be absorbed in the running quark mass. The same comment ap- 
plies to the contribution of "singlet" terms which are relevant for the axial part 
of the Z decay rate 0, |, 0, |. 

In principle one would of course prefer to control the mass dependence of the 
cross section for arbitrary m?/s. To date this has been achieved in first order 
in Us only. Higher orders in were calculated for the first term in the m? / s 
expansion, allowing for calculational techniques based on massless propagator 
type integrals. In this work this program will be pursued further and rn^ / 
terms will be calculated up to order a^. These terms are fairly unimportant for 
Z decays, however, they will become relevant in the low energy region. Examples 
are measurements at a i? meson factory just below or several GeV above the hb 
threshold, or alternatively around 5 GeV, an energy region that could be explored 
at the BEPC storage ring. 

The a1 calculation presented below is based on |10|. There the operator 
product expansion of two point correlators was studied up to two loops, including 
therefore terms of first order in a^. The expansion included power law suppressed 
terms up to operators of dimension four, which are induced by nonvanishing quark 
masses. Renormalisation group arguments, similar to those employed already in 
[§, |], will allow to deduce the a^m^ terms. The calculation is performed for vector 
and axial vector current correlators. The first one is of course relevant for electron 
positron annihilation into heavy quarks at arbitrary energies, the second one for 
Z decays into h quarks and for top production at a future linear collider. In order 
Us the terms can also be obtained from the Taylor expansion of the complete 
11,1^. This provides not only a (fairly trivial) cross check of the term. 



answer 



The comparison with the complete answer in first order may also indicate to 
which extent th? plus terms lead to a reliable estimate of the complete answer 
in order a1. 
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2 Mass terms of first order in ag 



QCD corrections to vector and axial current correlators in order as and for arbi- 
trary m^/s were derived in |]ll|. Compact formulae are given in [0. These are 



conventionally expressed in terms of the pole mass denoted by m in the following. 
It is straightforward to expand these results in m^/s and one obtains - in obvious 
notation - 
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The approximations to the correction function for the vector current correlator 
(including successively higher orders and without the factor ag/vr) are compared 
to the full result in Fig.|l]. For high energies, say for 2m},j \fs below 0.3, an excel- 
lent approximation is provided by the constant plus the vn? term. In the region of 
2m/ y/s above 0.3 the term becomes increasingly important. The inclusion of 
this term improves the agreement significantly and leads to an excellent approx- 
imation even up to 2m/y/s ~ 0.7 or 0.8. For the narrow region between 0.6 and 
0.8 the agreement is further improved through the term. Higher powers in m 
which can be calculated in a straightforward way, lead only to a modest further 
improvement in the region close to threshold, say above 0.8. This region with 
its behaviour governed by the Coulomb singularity requires therefore a different 
treatment. The same considerations apply for the axial current correlator (Fig.|I]). 

Inclusion of and terms hence allows to extend the QCD prediction 
from the high energy limit down to fairly low center of mass energy values, per- 
haps 4 to 5 GeV above the hh threshold, or from about 2 GeV above the charm 
threshold up to just below T(4S'). Similarly to the discussion in Q for the rn? 
terms, the logarithms accompanying the terms can also be absorbed through 



a redefinition of m in terms of the MS mass 13, 14 at scale s 
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and one obtains0(m = m(s)) 
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This resummation is possible for the second and fourth powers of m in first order 
Us and in fact for m? corrections to Ry and Ra in all orders of as- However, 
logarithmic terms persist in the rn^ corrections, starting from 0{a1), as discussed 
in the next section. Higher order terms in the mass expansion, starting from m^, 
evidently exhibit superficially "incurable" logs already in order as- In fact, these 
logs may be also summed up (see |l^, |10| and also [l^l)- Unhke the quadratic 



mass term, the resulting expression will contain some explicit dependence on the 
strong coupling constant as{mq) taken at the scale equal to the quark mass. The 
numerical effect of such a summation proves to be rather small unless the ratio 
s/m^ is chosen too large. 



3 Corrections of order ajm^ 

Motivated by the fact that the first few terms in the / s expansion provide 
already an excellent approximation to the complete answer in order a^, we now 
proceed to the evaluation of the three loop corrections. The calculation makes 
use of the properties of the operator product expansion as discussed in P, [10| . 
To allow for a comprehensive presentation, the relevant formulae are repeated 
whenever necessary. 

1. Operator Product Expansion (OPE) 

The calculation is based on the operator product expansion of the T-product of 
two vector or axial vector currents 

T,Aq, J) = ij T(J^(x) J+(0))e*^^da;, (6) 

with = U'-y^d (or = U'j^'y^d)- Here u and d are just two generically different 
quarks with masses rriu and rrid- Quarks which are not coupled to the external cur- 
rent will influence the result in order a^ through their coupling to the gluon field. 

•^The terms given below are in disagreement with those quoted in |]l5| . 
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Figure 1: Comparison between the complete 0{as) correction function (solid line) 
and approximations of increasing order (dashed lines) in for vector (upper 
Figure) and axial vector current (lower Figure) induced rates. 
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The result may be immediately transformed to the case of the electromagnetic 
current of a heavy, say, t (or 6 ) quark. 

The asymptotic behaviour of this (operator valued) function for = —q^ — > 
oo is given by an OPE of the form (Different powers of may be studied 
separately. Only terms proportional to Lorentz scalar operators of dimension 
four are kept to derive the corrections as eventually we are interested in the 
function flBD sandwiched between the vacuum states.) 



(7) 



As demonstrated in [0, ^ the following operators of dimension four may in 
general appear at the rhs. of (0), 



O5 = (V^^V)c^ 



r^ij 2 2 
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The (gauge non-invariant) operators O4 and O5 are non-physical; they appear 
since the gauge invariance of the Lagrangian is broken through gauge-fixing. Ma- 
trix elements of these operators vanish for physical and thus by definition gauge 
invariant states. Hence they do not contribute in the calculation below. They 
must, however, be taken into account in the calculation of the coefficient func- 
tions (CP's) multiplying the physical operators. Once these CP's are known, one 
may freely ignore the non-physical operators. 

In the following results were obtained for the (transversal) coefficient func- 
tions: 
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4L + ^(55-16C(3) + fL + 8L^ 



where L = In/x^/Q^, C(3) = 1.202. . ., / stands for the number of flavours; the 
upper (lower) sign refers to M7^(i (1*7^75^). 

From these results one may find for example the CF's for T^y{q, J) with = 
Irf^h or 67^756 through the substitutions d —>■ b, u ^ b and rrih, rriu — > w^fe. 

The starting point of the calculation was the nondiagonal ud current. Singlet 
contributions are therefore absent by construction. These would contribute in 
order a1 to the axial part of r(Z bb) and will not be considered in this work. 
Their relative importance will be discussed in section 4. 

It should be stressed that the generic OPE (0) determines the asymptotic 
behaviour of any Green function with an insertion of the T-product T^i,{q, J), in 
particular of the correlator 

U,,{q) = {T,,{q,J)). (12) 



In particular, as demonstrated in |]T0[ the expansion (|^ bears all the information 
about the large Q behaviour of the correlator (0) in perturbation theory, pro- 
vided the Vacuum Expectation Values (VEV's) of local operators are evaluated 
and renormalized according to the standard minimal subtraction prescription. 
(No normal ordering of the operator products!) The above results allow to find 
the absorptive part of the polarization operator Il{Q^) defined in the standard 
manner 

once the MS renormalized values for the vacuum expectations values (G^) and 
(qq) are known. 

However, this nice feature of the minimal subtraction scheme has its price: 
In schemes without normal ordering the renormalization properties of compos- 
ite operators are more involved. For instance, the operator miip^ipj mixes under 
renormalization group transformations with "operators" proportional to the unit 
"operator" times polynomials of fourth order in the quark masses. For our aims 
it will be necessary to know the anomalous dimension matrix of all physical op- 
erators in the list (|^). It has been shown in |T8|, |T0[ (for earlier related works see 



[p!9| , P0[ pi| , |22| ) that this matrix can be expressed, to all orders in perturbation 
theory, in terms of the QCD /^-function, the quark mass anomalous dimension 
7m («s) and the so-called vacuum anomalous dimension 70. The explicit expression 
for the anomalous dimension matrix will be presented below. 



6 



In the calculation of R{s) only the logarithmic parts of the CF's need to be 
taken into account (the constant terms have no absorptive parts at s > !). 
Three remarkable observations can be made: 

• In first order of ag only trivial operators (proportional to the unit operator 
times some combination of quark masses) contribute to R{s) since only 
these depend on L. Hence R{s) may be found directly from the coefficient 



function (11). This fact directly leads to the absence of the mass logarithms 
in R{s) in order ag as the latter may appear only from VEV's of non- 
trivial operators. This is in agreement with the expansion of the exact result 
discussed in section 2. Furthermore, the argument does not apply to terms 
of order a^m^, again in agreement with the expansion calculated in section 
2. 

• To find the correction of order a^m^/s^ to R{s) there is no need to compute 
three-loop contributions of order to the VEV's (G^) and (qq). Only 
already known two-loop terms of order (see below) are required. 

• The calculations of were performed at the two-loop level and thus their 
results as expressed by (^-(|ll|) do not include terms of order for the 
CF's Cg. Fortunately, the terms proportional to L may be inferred from 
the renormalization group invariance in analogy to [Q]. 

2. Renormalization Group (RG) Functions 
The following notation will be employed: 

• The effective couplant is defined as a{fi) = '^^^^^ and the number of fiavours 
is denoted by /. 

• The /3-function and the anomalous dimension of the mass 7m are defined 
through 

"/^ i>0 
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The anomalous dimension of the vacuum energy is0 
7o(a,m)=7o'(a)E"^' + 7o'(a)E 

i i^j 
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Operator mixing 



It has been found in [jlOl that the anomalous dimensions of the (physical) 
operators of dimension four contributing to (^ may be presented as follows^: 
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Ai = 47^06 



where the last equation is of course a direct consequence of the definition of the 
quark mass anomalous dimensions. 

4- RG constraints on the coefficient functions 

In our particular case of the polarization operator n((5^) the OPE (0) assumes 
the form 
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'''The leading term and the first correction in the expression for 7q were calculated in [|lO| . 
Recently one of us extended the calculation to include the correction of order a1 . 

^ The differences between ( pj| ) and eq.(3.4) of originate from a different overall sign in 
the definition of 7^ and the different normalization of the operator Qi. 
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Figure 2: Diagrams giving rise to nonzero vacuum expectation values of the com- 
posite operators ip'^' ci-nd G^^, in the lowest orders of perturbation theory. 



where 



(21) 



The condition of the RG invariance implies 



or, explicitly, 
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The structure of the last equation allows obviously to find the logarithmic terms 
of order a1 in Ce in terms of the function Ce and Ci which are to be known 
completely up to and including the order contributions and the coefficient 
functions C'g which should be known in the order aj.. An inspection of (|9|)-(|TlD 
immediately reveals that we indeed know these coefficient functions with the 
accuracy required. 

5. Calculation of VEV's of Oi and O2 

Diagrams contributing to VEV's of Oi and O2 in order are shown in Fig.0. 



The VEV of O2 was computed in ||2^, |T0[ and reads 

,2 / ,,2 
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The perturbative "gluon condensate" was computed in 



(eq. (B.la)): 
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We have confirmed the correctness of the result (^5[) by an independent calcula- 
tion. 



6. Calculation of the ajm corrections 

Now we are in a position to find the corrections of order alrn^ / to the both 
spectral functions Ry and Ra- Indeed, simply integrating the rhs of (|^) one gets 
the missing logarithmic terms of order a1 in Cq. 
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Note that the term proportional to J2i in the above result was independently 
found via a direct calculation in ||2f 



The next step is to use eqs. ( p4| , ^ , |26D to find the corrections of order a^m^/s^ 
to the spectral densities of the vector and axial vector correlators. The results 
reads (below we set for brevity the MS normalization scale /i = a/s and m„(s) = 
mrf(s) = m(s)) 
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Note that the sum over i includes also the quark coupled to the external current 
and with mass denoted by m. Hence in the case with one heavy quark u of mass 
m {d = u) one should set X^i =r = 1 and X]j =i- = 1- In the opposite case when 
one considers the correlator of light (massless) quarks the heavy quark appears 
only through its coupling to gluons. There one finds: 



Rv = Ra = +'^' 
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Numerically, eqs. 
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4 Discussion 



The Z decay rate is hardly affected by the contributions. The lowest order 
term evaluated with fn = 2.6 GeV leads to a relative suppression (enhancement) 
of about 5 X 10~^ for the vector (axial vector) current induced Z ^ bb rate. 
Terms of increasing order in become successively smaller. It is worth noting, 
however, that the corresponding series, evaluated in the onshell scheme, leads 
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to terms which are larger by about one order of magnitude and of oscillatory 
signs. The correction to r{Z qq) which starts in order is evidently even 
smaller. From these considerations it is clear that corrections to the Z decay 
rate are well under control — despite the still missing singlet piece — and that 
they can be neglected for all practical purposes. 

The situation is different in the low energy region, say several GeV above the 
charm or the bottom threshold. For definiteness the second case will be consid- 
ered and for simplicity all other masses will be put to zero^. The contributions 
to from terms are presented in Fig.^ as functions of 2m/ ^/s in the range 
from 0.05 to 1. As input parameters rripoie = 4.7GeV and Aj^g = 235MeV, cor- 
responding to as{m\) = -12 have been chosen. Corrections of higher orders are 
added successively. The prediction is fairly stable with increasing order in as a 
consequence of the fact that most large logarithms were absorbed in the running 
mass. The relative magnitude of the sequence of terms from the m? expansion 
is displayed in Fig.|. The curves for wP and are based on corrections up to 
third order in with the term starting at first order. The curve receives 
corrections from order zero to two. 

Of course, very close to threshold, say above 0.75 (corresponding to y/s below 
13 GeV) the approximation is expected to break down, as indicated already in 
Fig.|l|. Below the hh threshold, however, one may decouple the bottom quark and 
consider mass corrections from the charmed quark within the same formalism. 

Also Rqq where q denotes a massless quark coupled to the external current 
is affected by virtual or real heavy quark radiation. The m? terms have been 
calculated in and start from order a^: 

= -(-)' ^(15 -|/)(|-C(3)) (32) 
V 7r / s 3 3 

The terms of order a^m^ were given above. Both lead to corrections of (9(10~^), 
(evaluated at an energy y/s of 10 GeV) and can be neglected for all practical 
purposes. 

Summary: QCD corrections to the vector and axial vector current induced rates 
and cross sections of order a^m^ have been calculated. They have little effect 
on the Z decay rate. However, they are important for the analysis of hadron 
production in the low energy region. 

Acknowledgement We would like to thank A. Kwiatkowski and T. Teubner 
for help in preparing the figures. 

more detailed analysis of the phenomenological implications of these results for the low 
energy region will be presented in . 
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Figure 3: Contributions to BX from terms including successively higher orders 
in as (order a^/ a]/ corresponding to dotted/ dashed/ solid lines) as functions 
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